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A scaling law analysis of the world data on inclusive large-p^ hadron production in hadronic colli- 
sions is carried out. A significant deviation from leading-twist perturbative QCD predictions at next- 
to-leading order is reported. The observed discrepancy is largest at high values of x ± — 2p ± /y /r s. 
In contrast, the production of prompt photons and jets exhibits the scaling behavior which is close to 
the conformal limit, in agreement with the leading-twist expectation. These results bring evidence 
for a non-negligible contribution of higher-twist processes in large-p ± hadron production in hadronic 
collisions, where the hadron is produced directly in the hard subprocess rather than by gluon or 
quark jet fragmentation. Predictions for scaling exponents at RHIC and LHC are given, and it is 
suggested to trigger the isolated large-p^ hadron production to enhance higher-twist processes. 

PACS numbers: ll.15.Bt, 12.38.-t, 12.38.Qk, 13.85.Ni 



The production of a hadron at large transverse mo- where n spoc t a tor is the number of constituents of A, B, 
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mentum, p ± , in a hadronic collision is conventionally 
analyzed within the framework of perturbative QCD by 
convoluting the leading-twist (LT) 2 — > 2 hard subpro- 
cess cross sections with evolved structure and fragmen- 
tation functions. The most important discriminant of 
the twist of a perturbative QCD subprocess in a hard 
hadronic collision is the scaling of the inclusive invariant 
cross section [l|, Q , 

a ^^ E ^_ {AB ^CX) = ^^1, (1) 
d A p p™ 

at fixed x ± = 2p ± /y/s and center-of-mass (CM) angle 
z9. In the original parton model Q the power fall-off is 
simply n — 4 since the underlying 2 — > 2 subprocess am- 
plitude for point-like partons is scale invariant, and there 
is no dimensionful parameter as in a conformal theory. 
However, in general additional higher-twist (HT) contri- 
butions involving a larger number of elementary fields 
contributing to the hard subprocess, n ac ti V o > 4, are also 
expected. For example, the detected hadron C can be 
produced directly in the hard subprocess reaction as in 
an exclusive reaction. Such direct HT processes can give 
a significant contribution since there is no suppression 
from jet fragmentation at large momentum fraction car- 
ried by the hadron, z, and the trigger hadron is produced 
without any waste of energy. 

Apart from scaling violations due to the QCD run- 
ning coupling and the evolution of parton distributions 
functions (PDF) and fragmentation functions (FF), the 
invariant cross section of a given hard subprocess is ex- 
pected to scale quite generally as (neglecting spin correc- 
tions) d 



a inv (A B ^ C X) (x 
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and C not participating in the subprocess. From Eq. ©, 
it is clear that HT processes involving a large number 
of active fields will result into a p ± -exponent larger than 
the LT expectation {n > 4) , but will exhibit a slower fall- 
off with x ± from the smaller number of spectator fields. 
Therefore, at large x ± and not too large p ± , HT contri- 
butions to the cross section can become significant, lead- 
ing to an effective exponent higher than the LT expecta- 
tion. In Ref. [5| the cross sections of the HT subprocesses 
gq — > nq and qq —¥ ng, where the pion is produced di- 
rectly within the hard subprocess, have been calculated 
quantitatively in perturbative QCD and compared the 
results with the cross sections of the LT processes. This 
gives a contribution to Ea/d 3 p(AB 7rX) with nominal 
scaling n — 6 at fixed x ± and ?9 since n a ctive = 5 Q. 

In this Letter, the exponent n NLO of mid-rapidity 
particle production (t? = 7r/2) is computed in QCD 
at next-to- leading order (NLO) accuracy from Q, with 
CTEQ6.6 PDF and the de Florian-Sassot-Stratmann 
and Bourhis Fontannaz-Guillet FF into hadrons and 
photons 0, respectively. The x ± -dependence of n NLO 
at fixed p ± — 10 GeV is shown in Fig.[T]for pions, kaons, 
protons/antiprotons, and inclusive prompt photons. The 
hadron exponents increase slowly from n NLO ~ 5 at small 
values of x x (x ± = 1CP 2 ) up to n NLO ~ 6 at x ± = 0.5; 
there is very little dependence on the specific hadron 
species. The exponent extracted in the prompt pho- 
ton channel is below those of hadrons, by roughly one 
unit. The smaller photon exponent is understood from 
the (relative) absence of fragmentation processes and one 
less power in a s , leading to less scaling violation in this 
channel. Remarkably, n NLO is close to the conformal 
limit, n = 4, at the smallest values of x ± . 

In order to investigate possible HT dynamics in large- 
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71 data 


(n cxp ) 


(rz } 


E706 


M 


7T U 


31.6 


38.8 


2 - 


9 


10" 


l _ 


4 x 10" 


l 


25 


8.2 ±0.11 


6.1 ±0.09 


PHENIX/ISR 


[9,i0] 


7T° 


62.4 


22.4 


2 - 


7 


2 x 10" 


i _ 


2 x 10" 


i 


3 


7.5 ±0.19 


6.2 ±0.30 


PHENIX 




7T° 


62.4 


200 


2 - 


19 


7 x 10" 


2 _ 


2 x 10" 


i 


12 


6.7 ±0.05 


5.6 ±0.08 


UA1 


[15] 


h~ 


500 


900 


2 - 


9 


8 x 10" 


3 _ 


2 x 10 _ 


2 


18 


5.7 ±0.09 


5.2 ±0.04 


CDF 


[13] 


h ± 


630 


1800 


2 - 


9 


7 x 10" 


3 _ 


10~ 2 




5 


5.2 ±0.15 


5.0 ±0.07 




ri a~\ 
[14] 


tracks 


630 


1800 


2 


19 


7 x 10" 


3 _ 


2 x 10" 
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52 


5.7 ± 0.03 


5.0 ± 0.02 


CDF 


[16] 


7 


630 


1800 


11 - 


81 


3 x 10" 


•1 _ 


9 x lo- 


2 
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4.7 ±0.09 


4.3 ±0.01 


DO 




7 


630 


1800 


11 - 


107 


3 x 10" 


2 _ 


ur 1 




6 


4.5 ±0.12 


4.3 ±0.01 


CDF 


[19] 


jets 


546 


1800 


29 - 


190 


10" 


1 _ 


2 x 10" 


I 


9 


4.3 ±0.09 


4.6 ±0.01 


DO 


[20] 


jets 


630 


1800 


23 - 


376 


8 x 10" 


2 _ 


4 x 10" 


1 


23 


4.5 ±0.04 


4.6 ±0.01 



Table I: Data sets selected in the present Letter. The kinematical range P± in GeV), the mean (n exp ) extracted from each 
set composed of n data data points and the corresponding expectation in QCD at NLO, (n NLO ), are given. 
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Figure 1: x ± -dependence of n N for it (solid line), K 
(dotted), p/p (dashed) and 7 (dot-dashed), at p ± — 10 GeV. 

p ± hadron production, n exp has been systematically ex- 
tracted from measurements in p-p and p-p collisions, 
from fixed-target to collider experiments, and compared 
to LT QCD expectations. It is deduced from the com- 
parison of x ± -spectra at different CM energies, 

( J_ In (V5T/V55) ( ) 

which is equivalent to (JTJ) at fixed x ± . In order to re- 
duce systematic uncertainties, only experiments which 
measured x ± -spectra at two distinct CM energies are 
considered, except for the PHENIX results at y/s = 
62.4 GeV [9j compared to a fit of ISR measurements at 
yfs = 22 A GeV [lCf. The recent data analyzed in this 
Letter are summarized in Table HI The data sets include 
7T° measurements by the E706 at FNAL [ll[ and by the 
PHENIX collaboration at RHIC [l,[l3. At higher ener- 
gies, the measurements of charged hadrons (or charged 
tracks [11) inp-p collisions at y/1 = 630, 1800 GeV 
by CDPjia 0] and Vs = 500, 900 GeV by UA1 [H] 
are included in the analysis. Also considered are prompt 
photon and jet 0, H3 data obtained by CDF 



and DO at yfs = 546, 630, 1800 GeV. 

The hadron exponents plotted in Fig. [2] (left) exhibit 
a clear trend, with a significant rise of n exp as a function 
of x ± . Typical values of n cxp are n cxp ~ 5-6 at small 
x ± ~ 10~ 2 while PHENIX data point to a mean value 
n ox P = 6 7 ± o 05 at x ^ _ 1Q -i At higher values of 

x ± , the comparison of PHENIX with ISR data as well as 
the E706 measurements reveal an exponent even larger: 
n cx P _ 7 5 ± 0.19 ( x± = 0.2) and n cxp = 8.2 ± 0.11 
(x ± = 0.2-0.4), respectively. The E706 data clearly con- 
firm results reported long ago at the ISR, which are re- 
analyzed in a forthcoming paper [2lj ] . The results ob- 
tained in the photon and jet channels are strikingly differ- 
ent from what is observed for hadrons. Their exponents 
show almost no dependence on x ± , yet the data cover a 
wide complementary range: x ± — 0.04-0.1 for photons 
and x ± — 0.08-0.4 for jets. Importantly enough, the val- 
ues obtained lie only slightly above the conformal limit, 
n cxp ^ 4 g anc j exp ^ 4 mos ^. s ig n ifi can tly they are 

several units smaller than the hadron exponents taken at 
the same x ± (the p ± range being however different). 

In order to compare properly data and theory, NLO 
calculations have also been carried out within the same 
kinematical conditions as the experiments. The differ- 
ence between experimental and theoretical exponents, 
A(x ± ) = ri cxp — n NLO , is plotted in the right panel of 
Fig. [2] for hadrons and photons/jets. Note that the er- 
ror bars include both experimental as well as theoreti- 
cal errors, added in quadrature. The biggest theoreti- 
cal uncertainty comes from the variation of renormaliza- 
tion/factorization scales, for which all scales were varied 
from p^/2 to 2p ±1 as is common practice (the renormal- 
ization scale ambiguity can be removed using the meth- 
ods described in [22|). Fig. [5] (right) indicates that the 
hadronic exponents extracted experimentally prove sig- 
nificantly above the LT predictions. The discrepancy is 
moderate at small x ± , A(x ± ~ 10~ 2 ) ~ 0.5, but becomes 
increasingly larger at higher values of x ± : the PHENIX 
measurements at x ± ~ 10 _1 lead to A ~ 1 and the ex- 
ponent inferred from E706 data is two units above LT 
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Figure 2: Left: Values of n Gxp as a function of x ± for h /ir° (circles), 7 (squares) and jets (triangles). Right: A = n cxp — 1 
as a function of x ± , error bars include the experimental and the theoretical uncertainties added in quadrature (see text). 



expectations. In contrast, the scaling behaviors observed 
for photons and jets remarkably coincide, in excellent 
agreement with the NLO predictions. Part of the dis- 
crepancy between data and fixed-order calculations at 
large x ± ~ 1 could occur because of the appearance of 
large threshold logarithms. In (1 — x ± ), which should be 
resummed to all orders [231 ] ■ However, the discrepancy 
is also observed at small values oi x ± ~ 1CP 2 , where 
threshold effects are expected to be small. 

The most natural explanation for the hadron data is 
the presence of important HT contributions from pro- 
cesses in which the detected hadron appears in the hard 
subprocess. The dimension of the hadron distribution 
amplitude leads naturally larger exponents; see Eq. ©. 
In contrast, particles having no hadronic structure like 
isolated photons and jets are much less sensitive to such 
HT contributions and should behave closer to LT expec- 
tations, as observed. Another piece of evidence for HT 
effects is the larger exponents for protons than for pions 
observed at the ISR [2l| . According to Eq. © , the expo- 
nent of HT would be — 6 for pions (n ac tive = 5) and 



rip = 8 for protons (n a 
instead of n r , — n 



6), leading to n v 



at LT (see FigJT]). The experi- 
mental value obtained from the ISR [211 ] . n p — ~ I, 
thus reflects the mixture of LT and HT contributions to 
the total cross section. It has been noted [24| that the 
presence of color-transparent HT subprocesses such as 
uu — > pd can account for the anomalous features of pro- 
ton production seen in heavy ion collisions at RHIC [25j . 

In order to probe the HT effect explicitly, let us con- 
sider a 2-component model cross section with nominal 
power dependence 



niodcl 



(pp — > 7T X) OC 



A(x ± ) , B(x ± ) 



(4) 



corresponding to the LT (n activo = 4) and HT (n act i VG = 
5) processes, respectively. The actual p ± -exponents are 
modified by the running coupling and PDF and FF evo- 
lution. Assuming that the contributions to n NLO — 4 due 
to pQCD are the same for the LT and HT processes, Eq. 
(jl]) gives the effective exponent 



e{1 (x ± ,p ± ,B/A) 



diner 



model 



d\np 1 
2B/A 



+ n^(x ±lP J 



NLO 



B/A 



(x ± ,pj. (5) 



Note that n ot{ — > n 



NLO + 2 for B J A ->• 00. As shown in 
Fig. El (solid line), the LT pion exponent (evaluated at 
x x = 0.2) slowly decreases with p ± and reaches n = 4 
a,s p ± -> 00 because of asymptotic freedom. Eq. §5§ 
shows that rt off depends on the relative strength of HT 
corrections to the LT cross section, B/A. The value 
B/A ~ 50 GeV 2 is extracted from the data as shown by 
the dotted line in Fig. [3] However, a somewhat smaller 
estimate, B/A ~ 15 GeV 2 , is obtained when all scales 
are set to p ± /2 in the QCD calculation. We note that 
the HT rate for direct processes and therefore B/A are 
enhanced relative to fragmentation processes since the 
trigger hadron is produced without any waste of energy; 
thus the magnitude of the subprocess amplitude is max- 
imized since it is evaluated at the trigger p ± , and the 
initial momentum fractions x\ and xi are evaluated at 
small values where the PDF are largest. 

Finally, we discuss the phenomenological consequences 
of possible HT contributions to hadron production in p-p 
collisions at RHIC and LHC. In order to obtain qualita- 
tive predictions, the difference A between the experimen- 
tal and the NLO exponent has been fitted to the hadron 
data in Table [I] using a simple parametrization (with (p ± ) 
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Figure 3: p ± -dependence of n oft of pions at x ± — 0.2 in QCD 
at NLO (solid line) . The dotted line represents a fit based on 
a two-component model with B/A = 50 GeV 2 , see Eq. (0. 



the geometrical mean of the two experimental p^-bins) 



A™(x ± ,{pJ)=p (-logxJ Pl x 



2p 2 (l-xJ 



(p±) +p 2 (l-x x ) 



P3 



inspired by the 2-component model above described. 
As expected in QCD, A fit is vanishing in p ± —> oo 
limit at fixed x ± . This analytic form is somewhat ar- 
bitrary but flexible enough for making predictions be- 
yond the (x ± , p ± )-range probed in present experiments. 
The typical values of A expected at RHIC (taking 
= 200,500 GeV) and at the LHC {y/s = 7 TeV, 
compared to s/s = 1.8 TeV at Tevatron) are plotted as 
a function of x x in Fig. |U At RHIC, A fit is slightly be- 
low 1 at small x ± < 5.10 -2 but decreases towards zero 
at larger x ± (i.e. larger p ± ). At LHC, smaller devia- 
tions with NLO expectations are expected because of the 
large values of (p ± ) probed at high energy: A fit ~ 0.5 
below x ± = 5 x 10~ 3 (corresponding to p ± ~ 20 GeV 
at y/s — 7 TeV) and smaller above. From this, the ra- 
tios of x ± -spectra can be determined straightforwardly, 

R V^/V^ = (V^/v /i T) At,t+ " NL ° I where the NL0 expo- 
nents at RHIC (ti nlo ~ 5.3) and LHC (n NLO ~ 4.8) 
do not vary significantly in the considered x ± range. In 
order to enhance the HT contribution to hadron produc- 
tion, we suggest to trigger on isolated hadrons, i.e. with 
small hadronic background in their vicinity. The use of 
isolation cuts, usually applied for prompt photons, will 
strongly suppress LT processes. As a consequence, the 
scaling exponents of isolated hadrons are expected to be 
somewhat larger than in the inclusive channel. 

The evidence for higher-twist dynamics reported in 
this analysis supports the interpretation of heavy-ion col- 
lision measurements at RHIC, in which the dense QCD 
medium enhances HT contributions, and thus proton 




Figure 4: Predicted difference between the experimental and 
NLO scaling exponent at RHIC y/s = 200, 500 GeV and the 
LHC {\/~s — 7 TeV as compared to \fs = 1.8 TeV) based on 
a global fit of existing RHIC and Tevatron data 



production, by filtering LT processes due to partonic en- 
ergy loss 24J . Future RHIC and LHC measurements will 
provide further tests of the dynamics of large-p^ hadron 
production beyond leading twist. 
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